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Using three diflFerent Monte Carlo generators of high energy proton-proton collisions (HIJING, 
NEXUS, and PSM) we study the energy dependence of multiplicity distributions of charged particles 
' including the LHC energy range. Results are used for calculation of the information entropy, Renyi's 
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. dimensions and other multifractal characteristics of particle production. 
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nuclei is important issue of multiparticle dynamics 1] . Though MD contains information about particle correlations 

o 

in an integrated form it provides general and sensitive tool to probe the dynamics of the interaction. Analysis of the 

. . . 

I multiplicity data from pp and AA collisions will be thus important part of the physics programme at the Large Hadron 
Oh 

^ Collider (LHC) at CERN. In this new energy domain some basic questions concerning MD remain still open. Current 

^ ' theoretical understanding of the interplay between soft and semihard mechanisms of particle production is insufficient 

^ • to provide reliable estimates even of the elementary quantity characterizing probability distribution P„ of the n-particle 
_ _ I 

event - the average multiplicity of particles <n> produced in pp collisions. Different scenarios vary over a wide range 
of values 2]. Notwithstanding such uncertainty, large multiplicity of particles makes it feasible to study specific 
characteristics of MD which are relevant to fractal properties of the multiparticle dynamics. Fractality is usually 
connected with investigations of multiplicity in limited intervals of phase-space and demonstrated as a power-law in 



Energy dependence of multiplicity distribution (MD) of particles produced in high energy collisions of hadrons and 



resolution dependence of the multiplicity moments 



Such behaviour is characteristic for effects of nonstatistical 



fluctuations |5( and/or intermittency Besides fractality, there exists evidence in favour of multifractality 
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the structure of inelastic events. Some authors brought suggestions |9|| that the property could be expressed in simple 

n 

thermodynamic terms. Namely, it was shown 10] that constant specific heat, widely used in standard thermodynamics, 
reflects multifractal character of various stochastic systems in a reasonable approximation. In hadronic interactions, 
its constancy relative to the q-order number of the multifractal characteristics (the generalized dimensions) was 
indicated in Ref. 11 L Universality of the multifractal specific heat with respect to various hadron- nucleus reactions 



i i]. 



was demonstrated il^ by means of the method proposed by Takagi |13j |. Though exploiting similar variables and 
methods supports multifractal interpretation of data, conclusion about the methodology is far from being unique. In 
such situation, alternative approaches to the above aspects of multiparticlc production arc also needed. They can be 
helpful in better understanding of multifractality in high energy collisions and useful in extracting dynamical origin 
of its basic phcnomenological observation - the power law dependence of the corresponding measures with respect to 
the fractal resolution. 

In this paper we employ and further extend our earlier study of MD isj based on robust characteristics of MD 
such as Shannon's information entropy 



and its generalization, Renyi's order-q information entropy ^ 17 



(1) 



(2) 



for which Ii = S. They contain information on the multiplicity moments <n^>— ^ n^Pn m a non-trivial way. The 
outline of the paper is as follows. In Section 2 we show how regularities in energy dependence of the Shannon entropy 
provide independent constraint on the energy dependence of <n>. Section 3 is devoted to multifractal interpretation 
of observed regularities in terms of the Renyi's order-q entropy. We use three different Monte Carlo (MC) generators 
of high energy proton-proton collisions, HIJING 

Q, NEXUS QQ: 

and PSM |2l| to study energy dependence of 
MD up to the LHC energies. Results of MC simulations are reported in Section 4. Summary and discussion of the 
results is presented in Section 5. 



II. ENTROPY 



Entropy is important characteristic of systems with many degrees of freedom. It seems quite natural to use it in 
description of high energy multiparticlc production processes. In particular, entropy of MD is an effective variable 



characterizing inelastic collisions with many particles produced. The simple relation 

S -\n<n> ^- I il;{z)\nil;{z)dz (3) 



between the entropy S and the average multipUcity <n> is vahd with good accuracy for large enough <n>. The 
function •^'(z) is related to the probability P„ by the formula 

■0 (z = =<n> P„. (4) 

The KNO scaling Q of MD postulates energy independence of '>jj{z) what can be expressed according to Eq. Q as 
follows 

S — hi < n > ~ const. {^/s). (5) 

Last relation implies that all information about the energy dependence of MD is contained in its first moment <n>. 

II 

Such behaviour was indeed observed in hadron-hadron collisions up to the ISR energy range |23| . However, beyond this 
ran ge, t he KNO scaling of MD of charged hadrons produced in pp/pp collisions was shown to be significantly violated 
23. The observed scaling violation, i.e. break down of Eq. (|SJ) results via Eq. ||2Jl in non-trivial correlation 
between S and <7i>. Experimental situation concerning the difference — In <n> calculated for MD of charged 
particles in the full phase space is illustrated in Fig. 1. Here, due to the charge conservation, < n >—< nch/'2 > is 



half of the average multiplicity of charged particles. Recently published Tevatron data from pp collisions are also 
included in the figure. The values of the difference (|3J) show clear increase with energy and thus confirm violation of 
the KNO scaling observed at lower ISR energies. 

Let us note that increase of 5* — In < n > can not continue ad infinitum. It was shown ll^ that this difference is 
asymptotically bounded by unity from above. The extremal value of the r.h.s. of the relation (jJl is obtained from 
the principle of maximal entropy applied to all continuous functions ipi^)- If ttic functions I0J satisfy the conditions 

/ tp{z)dz = / ztp{z)dz = 1 (6) 
Jo Jo 

valid for all normalized MD with given <n>, the extremum is reached for the KNO function tpiz) — exp(— z) which 
corresponds to the geometrical distribution. The difference 

S - In <n>< (1+ <n>)ln(l + = 1 + 0{<n>-^) (7) 

<n> 

is thus bounded from above by the expression which is calculated for the geometrical MD. The bounding values are 
represented by the curve in Fig. 1 . Last relation shows that energy dependence of the information entropy can be used 
as an independent constraint on energy dependence of the average multiplicity. 



p. 



Energy dependence of the entropy S in high energy coUisions was first studied in Ref. '141. Using MD of charged 
secondaries produced in pp and pp collisions in the energy range y/s < 900 GeV, monotonous increase of S with 
center-of-mass energy ^/s was found. For ^/s > 20 GeV, the linearity 

S = D^Yrn (8) 

with the maximum rapidity Ym = ln(y^/TO^) of the hadrons produced is valid. Here ttIt^ is the pion mass and Di is 
an energy independent constant. Recently published Tevatron data [2^ extend the validity of these findings up to 
= 1.8 TeV. The experimental situation is shown in Fig. 2 where the published errors of MD have been taken into 
account. Value of Di is within ±2% error band consistent with the predicted one Di = S/Ym = 0.417 ± 0.009 ij 
shown by the full/dashed lines in Fig. 2. On closer inspection, however, one can see that Di calculated form E735 
data is systematically above the UA5 values. It is connected with higher tails of the E735 MD as compared to the 
UA5 data. The discrepancy might be due to the sizable systematic uncertainty in both E735 and UA5 since in both 
experiments the full phase space MD were obtained by computer simulation from data measured in a restricted range 
of rapidity. 

As the bounding value in lO tends to unity, the observed monotonous increase of the entropy S = Di \n[^/'m-^), 
if valid at higher energies, will also govern the energy dependence of <n>: 

<n>K exp(S') = iV^/m^)^' . (9) 



Such asymptotic power law behaviour of the average multiplicity should be contrasted with other approaches to 
multiparticlc production. In particular in Ref. Q], <n> was predicted to increase as a second order polynomial in 
Ins. Difference between these two predictions will be substantial at the top energy of LHC, because according to Eq. 
© <n>f» 110 at ^/s = 14 TeV, while according to the parametrization used in Ref. j2|| <n>~ 70 at this energy. 

III. MULTIFRACTALITY 

Fractal geometry is nowadays widely used in many branches of physics 17 1 . In astronomy analysis of galaxy and 
cluster distributions it has led to surprising result that galaxy correlations up to 150h~^Mpc are scale invariant and 
not homogenous 27 1 . 
of stochastic systems 



In multiparticlc dynamics methods introduced originally for description of the fractal properties 

y. 



28l| are used extensively 



In particular, study of MD in small rapidity bins using the scaled 



factorial moments Fq revealed typical linear behaviour 

\nFq = ^TqlnS (10) 



in terms of the bin resolution S. Such behaviour was interpreted as (multi)fractal property of particle production. 
Another approach to the multifractality is directly connected to the Renyi's order-q information entropy Iq 1^)1. 
Energy dependence of Iq in high energy collisions was first studied in Ref . ^ . Using MD of charged secondaries 
produced in pp and pp collisions in the energy range ^/s < 900 GeV it was found that increase of Iq with center-of-mass 
energy is similar for various values of q. For ^/s > 20 GeV, the observed asymptotic linearity 

Ig = DqYra (11) 

with maximum rapidity of the hadrons produced, = lii{y/s / 171^^) , generalize Eq. Let us briefly summarize 

multifractal interpretation of Eq. In contrast to standard intermittency and multifractal analysis we consider 

fractal resolution S related to the total energy -^/s available and not to the phase-space bining. During hadron-hadron 
interaction with many particles produced, the energy dissipates into iV — ^/s/mT^ discrete sites each of the size TOtt. 
The site labeled by n is occupied with the probability P„ . Since most of the sites are unoccupied the overlay of many 
inelastic events can be visualized as a fractal with overall extent ^/s characterized by a local mass distribution function. 
The sufficient condition to produce such self-similar (hierarchical) structure is that the probability P„ exhibits some 
type of scale invariant behaviour. In particular, if at sufficiently high energies, i.e. at high enough resolution 5 — 1/N , 
the Pn acquires a power law dependence on the resolution 5, the quantity "^{PnY '^i^ scale with 5 like 



^(p„)«^r(i-9)^. (12) 
[7 



If D > is independent of 5, one usually speaks about fractality [I3 of the distribution P„. The multifractals 
[17 , 2^ generalize the notion of fractals for D = Dg depending on q. Spectrum of generalized dimensions Dq which 



has the following meaning. Dq corresponds to the capacity (box 
1 2 characterizes scaling of the information 

[33 



for multifractals is a decreasing function of q 

dimension) of the support of the measure P„, information dimension D 

entropy and Dq^s for integer q > 2 can be related to the scaling behaviour of g-point correlation integrals 
The observed approximate independence of the generalized dimensions 



Dq = -Iq/llld (13) 

n 

on the energy (resolution S) as well as their decrease with increasing q show 15] that for y/s > 20 GeV the full phase 
space MDs of charged particles from non-single-diffractive hadron-hadron collisions are indeed multifractal. Let us 
note that multifractality besides predicting Dq to be decreasing functions of q does not, in general, provide any further 
information about the g-dependence of the spectrum of the generalized dimensions Dq. In particular knowledge of 



say Di and D2 is insufficient to predict scaling behaviour of the higher q correlation integrals. It is thus gratifying to 
find out that this could be, at least in principle, achieved within interpretation of multifractality in thermodynamical 
terms [10|. Latter is based on analogy between l.h.s. of Eq. l(T^ and partition function 



Z(g)^^(P„)'' (14) 



with q playing the role of inverse temperature q = T^^ and V = — In 5 representing volume. The thermodynamic 
limit of infinite volume y ^ cxi is then equivalent to the limit of increasing resolution S 0. In the constant specific 
heat approximation the q-dependence of the generalized dimensions Dq acquires particularly simple form 



im s 
3 



„ . N Ino , , 

Dg^{a-c) + cj-^. (15) 

The coefficient c represents multifractal specific heat and a = Di . Regular behaviour of this type is expected to occur 
for multifractals for which, in classical analogy with specific heat of gases and solids, the multifractal specific heat c 
is independent of temperature [2^ in a wide range of q. 

We have examined validity of the approximation given by Eq. (|15|l for Dg defined by Eqs. Hll|) and l(2Jl. In Fig. 3a 
we present q-dependence of generalized dimensions calculated from the Tevatron data at y/s = 300, 546, 1000, and 
1800 GeV. One can see from the figure that the values of Dq reveal indeed linear increase as a function of \n{q)/{q— 1). 
This behaviour makes it possible to define the slope parameter c in the region q > I. Similar Dq dependencies for 



033 



are shown in Fig. 3b. The dashed line coincides 



data from CERN ISR and SppS Collider experiments 
with the full line in Fig. 3a indicating position of Dq values calculated from E735 data. Both data sets obtained by 
Tevatron and CERN experiments reveal approximately the same slope while their intercepts are mutually shifted. 
The shift is due to larger values of Di for E735 data, as already shown in Fig. 2. As pointed in the previous section, 
the discrepancy in intercepts is connected with systematically larger high multiplicity tails of the data from Tevatron 
when compared to the data from the CERN SppS Collider. This might be connected with the mutual systematic 
uncertainties of the experimental procedures when extending measured data into the full phase space region. 

Fitting the slope parameters for g > 1 at each separate energy, we have determined the values of multifractal 
specific heat c. The results are presented in the lower part of Fig. 2. The errors in determination of c were calculated 
from error bars of MD quoted in literature. They represent ^ 10 — 15% of the established values. For y/s > 20 GeV, 
the multifractal specific heat is within the estimated errors approximately energy independent quantity and achieves 
the value c « 0.08. Note that this number practically coincides with the slope obtained from the electron-positron 



multiplicity data 



3] 



what is indicated by the full line in Fig. 3b. Since our study concerns the full phase space and 



it is performed in a different sense than the usual intermittency analysis, we obtain smaller value of c in comparison 
with the specific heat (c ~ 0.26) determined from multifractal properties of the factorial moments While the 

multifractal specific heat reported in Ref. \l\ reveals some kind of universality with respect to various interactions, 



energy dependence of Dq obtained with the same method Il3| seems to be significant. Contrary to this, our method 
gives smaller values of the generalized dimensions Dq which are approximately energy independent for > 20 GeV. 

Energy independence of the multifractal specific heat confirms that, in addition to the information dimension Di, 
there appears to be yet another parameter c which could be used as universal characteristic of particle production in 
hadron-hadron interactions at high energies. Knowing Di and c, all other Dq can be thus deduced from Eq. 1151 

IV. MONTE CARLO SIMULATIONS 

We have exploited three different Monte Carlo models in our investigation of the proton-proton interactions beyond 
up to date accessible energies. 

n 

1. The HIJING 18] Monte Carlo code is based on QCD-inspired models for multiple jets production. It allows to 
study jets and the associated particle multiplicities. The model includes minijet production, soft excitation, nuclear 
shadowing of parton distribution functions and jet interaction in the dense nuclear matter. 



2. The NEXUS 



Q 2^ 



Monte Carlo code is based upon the hypothesis that the behaviour of the high energy 
interactions is universal. Basic building blocks of hadron-hadron or nucleus- nucleus scattering are parton ladders 
coupled softly to the nucleons. It relays on a consistent multiple scattering approach in the sense that most of the 
dynamics follows from a formula for the total cross section expressed in terms of cut diagrams. 

3. The Parton String Model (PSM) ^21j includes in its initial stage both soft and semihard components which 
lead to the formation of color strings. Collectivity is taken into account considering the possibility of strings in color 
representations higher than triplet or anti-triplet by means of string fusion. String breaking leads to the production 
of secondaries. 

We have analyzed y^-dependence of the average multiplicity <nch>, the entropy S, the difference S —\n{<nch> /2), 
and the multifractal specific heat c for the charged particles in the energy interval 25 GeV< y/s < 14 TeV. The 
quantities have been obtained for 14 points from this interval using 10000 simulated events for each one. The results 
of the energy dependence of <nch> are shown in Fig. 4. Though all three models give similar values in the range 
y/s < 1 TeV, their predictions significantly differ for y/s >1 TeV. The average multiplicity simulated by HIJING and 
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PSM codes is well approximated by the power dependence 

< rich >= as'' (16) 

with b ~ 0.018 for HIJING and PSM. The NEXUS predictions faU weU bellow and can be fitted much better by the 
second order polynomial in logarithm ^/s: 

< rich >= cio + ai In ^ + a2 In^ y/s. (17) 

The above parametrization corresponds to the prediction of the average multiplicity given in Ref. 

Let us now compare both parametrizations from the point of view of the observed regularity in the information 
dimension Di. The energy dependence of Di = S/Y^ is displayed in Fig. 5. For HIJING and PSM the behaviour 
is consistent with the experimentally observed ratio S/Ym = 0.417 ± 0.009 extrapolated towards super-high energies. 
This supports the conjecture that information entropy per unit rapidity should stay constant in the full phase space 
up to the LHC energy region and even beyond. In such case, energy dependence of the average multiplicity should 
achieve the asymptotic power law behaviour The HIJING and PSM simulations of the average multiplicity Hlt)|) 
confirm this tendency. On the other hand, results of simulations with the NEXUS code show slow but continuous 
decrease of Di = S/Y„i with the energy. The decrease of Di is connected with non-power like, namely the logarithmic 
energy dependence H17|) of the average multiplicity. Using the general limit iQ on MD one can say: either the entropy 
S must slow down and violate the entropy scaling ijSJl as indicated by the NEXUS model, or, providing the entropy 
scaling stays valid, the average multiplicity must grow faster, similar as in the HIJING and PSM models, and violate 
thus the parametrization 1)17(1 of Ref. 

The above statement relies in details on the energy dependence of the difference 5 — ln(< rich > /2) which is depicted 
in Fig. 6. Though for all energies considered the difference does not reach the bounding value Q, its saturation in 
the LHC energy range is almost complete for the PSM model. The NEXUS prediction falls well bellow at the LHC 
energy and is closer to the KNO-scaling behaviour given by Eq. (|SJ). In the case of the HIJING we have checked the 
sensitivity of S* — ln(< rich > /2) to the interactions dynamics, changing in particular number of the produced jets 
Njcts- There are points of regime changes clearly visible in Fig. 6. They correspond to various maximal number of 
produced jets, which had been set to Njcts =0,1, and 5 for yG < 0.1 TeV, 0.1 TeV< y/s < 1 TeV, and y/s > 1 TeV, 
respectively. Character of these changes demonstrates that the quantity S — lii{<nch> /2) is sensitive to the dynamic 
of the interaction, in particular, to the number of produced jets. Increase of Njcts results in larger number of parton- 
parton collisions which play role in the KNO scaling violation 26]. The violation seems however not to destroy the 
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self-similarity of parton dynamics as manifested by the regularity in behaviour of the information entropy of MD 
including the data from Tevatron (Fig. 2). 

Multifractal character of particle production consistent with approximate energy independence of the multifractal 
specific heat c is predicted by the PSM, HIJING, and NEXUS models in the LHC energy region and beyond. We 
demonstrate this in Fig. 7. The tendency towards a constant value of c confirmed by the simulation results is similar 
with the experimental situation shown in Fig. 2. It will be interesting to see whether this will really hold true for the 
forthcoming data from LHC. 



V. SUMMARY AND DISCUSSION 

Having enumerated experimentally observed features of MD in multi-hadron production such as regularities in 
the information entropy and generalized dimensions, we have studied the multifractality of inelastic hadron-hadron 
collisions at high energies. We have applied the multifractal specific heat approach to characterize this property in 
a quantitative way. The generalized dimensions were calculated using the resolution S — raT^j ^ defined as minimal 
fraction of the center-of-mass energy which can be carried away by the outgoing hadrons. The experimental data 
on MD are consistent with the constant specific heat approximation c ^ 0.08 with respect to both, the multifractal 
temperature T = ^ and the energy ^/s. 

We have extended our study of MD in pp collisions up to LHC energies exploiting three different Monte Carlo 
generators of high energy hadron-hadron and nucleus-nucleus collisions (HIJING, NEXUS, and PSM). The HIJING 
and PSM predictions confirm approximate energy independence of the entropy dimension Di = S/Ym = 0.417±0.009 
\l\ in the full phase-space up to the LHC energy and even beyond. These models give Slhc ~ 4.65 at -^5=14 TeV 
and a power-like energy behaviour of the average multiplicity. Simulations with NEXUS show monotonous decrease 
of S/Y„i which corresponds to a slower (logarithmic) increase of the average multiplicity with energy. In this view, 
HIJING and PSM, unlike to NEXUS, prefer self-similar character of multi-hadron production. All three models 
predict approximately constant value of the multifractal specific heat and thus suggest multifractal character of 
proton-proton interactions in the considered energy region. The information about self-similarity and multifractal 
structure is contained in MDs and can be quantified by the behaviour of both the information entropy and the 
multifractal specific heat in a model independent way. 
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FIG. 1: The energy dependence of the difference S — ln(<nc/i > /2) for charged particles (experimental data). The full circles 
correspond to data from E735 experiment. 
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FIG. 2: The energy dependence of the entropy dimension Di and the multifractal specific heat c for charged particles. The 
triangles correspond to data form E735 experiment. 
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FIG. 3: The generalized dimensions Da as function of ln(q)/(q-l) for charged particles. Data are taken (a) from E735 



experiment |2q | and (b) from Refs. 
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FIG. 4: The energy dependence of the average multiphcity simulated by MC for charged particles. 
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FIG. 5: The energy dependence of Di = S/Ym of charged particles. The points correspond to the MC simulations. 
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FIG. 6: The energy dependence of the difference S — ln(< rich > /2) of charged particles. The points correspond to the MC 
simulations. 
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FIG. 7: The multifractal specific heat c calculated from MC simulations. 



